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Abstract 

The quantum yield of semiconducting single walled carbon nanotubes is limited by strong 
diffusion driven end/defect quenching and rapid Auger-like exciton-exciton annihilation (EEA) 
processes. While these two processes have been extensively studied independently, we show 
that both effects must be considered to properly extract intrinsic nanotube properties. To 
achieve this result, we have developed a convenient analytic model of exciton dynamics (ver- 
ified by comparison to Monte Carlo simulations) that shows the intrinsic optical quantum ef- 
ficiency is substantially higher than previously reported from observations in typical samples. 
Diffusion limited EEA is modeled using a one dimensional line segment revealing enhanced 
annihilation at low exciton populations. Lastly, the analytic model is demonstrated using non- 
linear EEA effects to extract intrinsic absorption cross sections, lifetimes, and diffusion lengths 
for SWCNTs. 
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Introduction 



Many interesting photo-physical applications utilizing single walled carbon nanotubes 
(SWCNTs) are emerging that take advantage of two dominant non-radiative decay pathways, 
end/defect quenching and rapid Auger like exciton-exciton annihilation (EEA). 1 At low exci- 
ton densities, the stochastic fluctuations of the optical emission intensity due to adsorption and 
desorption of fluorescence quenchers has led to single molecule counting arrays and biological 
sensors.- - - At high exciton densities, the rapid relaxation of excitons has been useful in creat- 
ing saturable absorbers for passive mode-locking 6 " 9 and may give rise to photon anti -bunching,— 
showing promise for single photon sources with selectable wavelengths from optical frequencies 
into the telecom band. 

Following the work of Cognet and Weisemann, first demonstrating stochastic quenching of ex- 
citons in SWCNTs, 1 1 much effort has gone into modeling the length and defect dependent quantum 
efficiency and exciton lifetimes.—"— These models have been developed ignoring EEA since it was 
expected that at low excitation fluences non-radiative decay would be dominated by end/defect 
quenching. 12 However, the fluence at which EEA plays a role is surprisingly low, as has been 
observed by nonlinear photoluminescence (PL) effects for both continuous wave (CW) 18,19 and 
pulsed^" 2 ! laser excitation, putting into question whether EEA can be ignored for low to moderate 
fluences. 

It is clear that as the excitation fluence is increased, EEA dominates the exciton decay dynamics 
and substantially reduces the optical quantum efficiency. The rapid injection of many excitons un- 
der pulsed excitation leads to a complete saturation of the SWCNT photoluminescence- 2 ^ 2 ! while 
even under CW excitation a strong nonlinearity has been observed at surprisingly low fluences— ^ 
and was proposed to be an effect of long lived optical quenchers. 19 For both pulsed and CW excita- 
tion, the nonlinear fluence dependence of the photoluminescence has been modeled as a diffusion 
limited process.— £222. However, a full microscopic model of exciton dynamics which properly ac- 
counts for the spatial and density dependence of exciton lifetime is still absent. In fact, the existing 
analytic diffusion limited EEA model 23 approximates a SWCNT as a ring which leads to an incor- 
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rect exciton pair distribution function, overestimating the exciton separation and underestimating 
EEA at low exciton levels. In this work we consider three competing decay processes, linear decay 
(Pr), decay due to exciton-exciton annihilation (Peea), an d decay due to end quenching (Pend)- 
This microscopic model can be used to accurately determine intrinsic SWCNT properties such as 
the optical absorption, quantum efficiency, exciton lifetime, and diffusion length for any SWCNT 
sample. 

Results and discussion 

Exciton End Quenching 

While the study of exciton diffusion limited quenching at SWCNT ends and defect sites has been 
explored extensively, 13 i 15 ~ i& it has not been developed considering many body exciton-exciton 
annihilation and does not properly account for the spatial dependence of the effective exciton 
lifetime. For this reason, the following diffusion-end quenching formalism is developed which can 
be easily incorporated into a diffusion limited EEA model. An exciton can be treated diffusively 
and not ballistically as long as, L >> y/2DxAt, where L is the nanotube length, Dx is the diffusion 
constant, and At is the scattering time and is approximately equal to the dephasing time T2.— In 
this case, the diffusion of an exciton in one dimension obeys the diffusion equation 8 t n(x,t) = 
Dx5^n(x,t). If the SWCNT ends quench with a probability of 1, then the probability that an 
exciton created at a position x on a line segment [0,L] will survive after time t is equal to the 
probability that it has not passed through the positions a distance x or L — x away. From the first 
passage time of a Brownian motion, 24 the probability that an exciton, starting at position 0, has not 
passed a position a distance x away after time t is, 




(1) 



where erf(z) = (2/y/n) J e ' dt is the error function. 
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For a SWCNT of length L with only the two ends as quenching sites, such as suspended from 
a trench, the average probability that an exciton survives after time t is, 

L 

J g(x)P(x,t)P(L-x,t)dx 

Pend (L,t) = ^ - L (2) 

Jg(x)dx 
o 

where g(x) is the spatial weighting function determined by the excitation source. In the case of 
homogeneous excitation, g(x) = 1 resulting in, 



END 



erf 




32? 
%x D 



V2 



erf 




erf 




e "sr 



(3) 



where Xq = L 2 /Dx is the only physical parameter and is essentially the average "diffusional 
time". 23 In the limit that t « Xd Eq. (3) simplifies to Pend^dJ) ~ 1 — A^Jt/nxo but in general 



calculations are not done in this limit since we are interested in both short and long time dynamics. 
The exciton population decay due to end effects [Eq. (3)] is shown in Figure 1(a) for several values 
of Xr>- It is useful to note that the exciton population decay to the edges without intrinsic decay is 
on time scales comparable to those measured in PL lifetime measurements.— In general, because 
of the spatial dependence of the exciton population decay, this is not equivalent to a Poissonian 
decay with an effective decay rate T e ff, as has commonly been modeled. 

In addition to non-radiative end quenching excitons also undergo linear radiative and non- 
radiative decay with the survival probability Pr(t) = e~ Tt . Here r = Tr + T^r is the linear decay 
rate from radiative (Tr) and non-radiative (T^r) linear decay. End quenching is taken to be in- 
dependent of linear decay so the survival probabilities can be multiplied and the exciton decay 
dynamics follow, P[xo,t) = Pend(^dJ) x Pr(t). To determine the quantum efficiency (t\qe) of a 
system for which the exciton population decay follows an arbitrary function, we note that the prob- 
ability of radiative decay in an infinitesimal time dt is YRdt. It follows that the number of excitons 
that radiatively decay at time t is simply the probability of radiative decay times the total number 
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of excitons at time t, N ra d(t) = N ex {t)T R dt where N ex (t) = NoP(t) starting with the injected exci- 

oo 

ton number No. The total number of excitons that radiatively decay is then N ra d = NqT r J P(t)dt, 

o 

allowing the quantum efficiency can be written as, 



i Q E = ^- = r R J p(t)dt (4) 



for any arbitrary decay function P(t). In the simple case where end effects are negligible 
(for example L — > oo) Pend(°°J) — > 1 and the exciton population decay (Pr x Pend) simplifies 
to a standard mono-exponential decay, 26 P(t) = e~ Tt , and the solution to Eq. (4) is the intrinsic 
quantum efficiency rjo = T R /T. 

Typical SWCNT lengths range from ~ 100 nm to > 50 /im— >2Z and diffusivities range from ~ 
0.1 - 100 cm 2 /s.— Within this typical range, there is over an order of magnitude decrease in 
the optical quantum efficiency compared to the intrinsic quantum efficiency due to end quenching 
as can be seen in Figure 1(b). As verification, the analytic model is compared to standard Monte 
Carlo simulations [see Methods], which are processor and time intensive and are of limited value 
in identifying the importance of the underlying microscopic processes. The effect of a SWCNT 
ensemble with random defect sites can easily be calculated by averaging over the permutations of 
smaller effective lengths. 




Time (ps) t d (ns) SWCNT Length (nm) 

Figure 1: (a) Exciton survival probability (Pend) due to end quenching [Eq. (3)] for several values 
of the average diffusional time Tq. (b) Quantum efficiency [Eq. (4)] as a function of the average 
diffusional time (To) with T = 3.3 ns -1 . (c) Photoluminescence from "pristine" SWCNTs (r = 
3.3 ns~ ! and Dx = 44 cm 2 /s) of various lengths excited with a Gaussian beam (o = 520 nm) and 
[INSET] compared to experimental data collected by Moritsubo et al— 
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As an example of the utility of this formalism, the photoluminescence dependence on tube 
length can be calculated since PL T}qe, shown in Figure 1(c). This can be compared to exper- 
imental results collected for variable tube lengths (such as done by Moritsubo et al.— ), shown 
in Figure 1(c) [INSET], and can be used to extract intrinsic properties. Experimental studies 
use Gaussian beam excitation (a = 520 nm) which are included using Eq. (2) with g(x) = 
exp( — (x — L/2) 2 /o 2 ). The measured 8 meV linewidth corresponds to a room temperature thermal 
diffusivity of Dx = 44 cm 2 /*.— It should be noted that in general, intrinsic lifetimes (i.e.,T = 1 /T) 
is difficult to directly obtain from experimental time-resolved population measurements, since such 
studies measure effective lifetimes which include quenching at ends. Moritsubo et al. use 85 ps 
effective lifetimes based on single SWCNT femtosecond excitation correlation (FEC) measure- 
ments.— Analysis of the original FEC data using the model outlined in this manuscript, extracts 
an intrinsic effective lifetime around four times longer than the effective value (t = l/T = 300 
ps). With one additional fit parameter (for y-scaling), the analytic model prediction shows good 
agreement with experimental observations (Figure 1(c) [INSET]). Longer SWCNT lengths would 
further constrain parameter fitting. 

The intrinsic exciton diffusion length \/Dxi is found here to be ~ 1.2 fim. For finite length 
SWCNTs, the observed diffusion length is shorter since the effective lifetime is shortened by end 
effects. Though this analysis allows extraction of certain SWCNT parameters such as diffusion 
length, results are still linearly dependent on exciton number and therefore cannot decouple the 
PL action cross section (PLAC=T]g£ x o a b) to extract the quantum efficiency (t\qe) and the optical 
absorption cross section (o a b)- To do this, experimental results at higher fluences involving exciton- 
exciton annihilation are required. 

Exciton-Exciton Annihilation 

As the excitation fluence is increased, for both CW— ^ and pulsed 20,21 excitation, EEA substan- 
tially reduces the optical quantum efficiency. Even with average exciton densities less than one 
exciton injected for every diffusion length, we show that diffusion mediated annihilation can oc- 
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cur. 

Previous work on exciton diffusion-reaction modeling 23 in carbon nanotubes used the geome- 
try of a one dimensional ring 30 which is a good approximation for infinitely long SWCNTs and/or 
large exciton densities, but is not in general correct. In many SWCNT samples the ratio of diffu- 
sion length to tube length can be quite high 13 -i£ri£ and PL saturation has been observed even for 
quite low exciton densities making this regime important to model correctly. 21 The assertion that 
the pair distribution function for the ring model is not appropriate can be easily understood from 
a two exciton example; in a ring geometry all possible initial pair separation lengths are equally 
likely 30 while for a line segment there are many ways to place an exciton pair separated by a few 
nanometers but there is only one way to place two excitons the length of the nanotube apart. Since, 
the excitons on average will be farther apart for the ring pair distribution function, application of 
this model will underestimate the role of EEA. To model the dynamics on a finite line segment, the 
pair distribution function for N particles is,— 



f(L x ,N)dL x = - 

^^(L-L x f- l dL X: R^0 

where L is the SWCNT length, R is the exciton size, L x is the nearest neighbor exciton separa- 
tion, and 6 = NR/L. In this work, we treat the exciton as a point particle (R — >■ 0) which simplifies 
Eq. (5). As a simple example of the difference between ring and line models, the pair distribution 
function for a line has a mean separation of L/(N+ 1) and not L/N as expected for a ring.— In 
the thermodynamic limit the pair distribution functions converge, but for low exciton densities, 
the focus of this work, the geometrical distinction is significant and the annihilation dynamics are 
different. 

Before expanding into the N exciton case, we first consider the simple case of two excitons with 
an initial separation L x independently undergoing one dimensional diffusion. In this model two 
excitons are taken to annihilate if they spatially overlap in their lifetime, resulting in one exciton 
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being destroyed and the other promoting to a higher energy state from which it quickly relaxes back 
into a thermal En state. Similar to a single exciton undergoing diffusion to a stationary quenching 
site (Eq. (1)), the probability that the separation (initially Lx) of two excitons undergoing relative 
diffusion has not reached after time t can be written as, erf (Lx / (2\/2Dxi)) . This is equivalent 
to a single exciton diffusing with diffusion constant 2Dx- Therefore the average probability that 
an exciton pair has not undergone EEA after time t is the mean of survival probabilities weighted 
with the pair distribution function (Eq. (5)): 




For the case of N = 2, this simplifies to, 

where again, the only physical parameter is To, the average diffusional time. This contrasts 
with the ring model which for two particle EEA predicts, P^ A (2, XdJ) = er f In 
Figure 2, the results of the two models (for To = 5.7 ns) are compared to Monte Carlo simulations 
for exciton dynamics in a finite length SWCNT. It can be seen that for the two exciton case the 
exciton pair survivability for typical exciton lifetime (~ 50-300 ps) is reduced by ~ 30% compared 
to predictions from previous work. 

Using the exciton pair survivability (Eq. (6)), the total number of excitons that survive EEA 
after time t, can be written as (l + (N — \ ) P^ A (N ,XD,t)), which has been verified by Monte 
Carlo simulations. This assumes one exciton of the pair survives EEA. The solution for this is the 
generalized hypergeometric function, 
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Figure 2: Analytic solutions for the EEA pair survival probability for a ring u and line segment for 
a 5 jUm SWCNT with Dx =44 cm 2 /s compared to Monte Carlo simulations. 
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(8) 



where r(jc) is the Gamma function. From this it follows that, given Nq initial excitons, the 
probability an exciton survives EEA is simply Peea(No, Tz>, t) = N ex (No, To, t) /No. In the case of a 
single injected exciton, the survival probability Peea remains 1 for all times t as expected. 

Combined Exciton Decay Dynamics 

Finally, to model the decay dynamics of the entire system including the three decay channels 
(PfiPeeAiPend), the exciton survival probabilities of the three decay channels are treated as inde- 
pendent so the total survival probability of the exciton population is simply, 



PtOt{N , H),r,0 = Pt(T,t) X PEND(lD,t) X P EE a(N , % D ,t) 



(9) 
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This model ignores the interplay between Peea and both Pend an d Pr, which slightly overes- 
timates EEA and consequently underestimates the quantum efficiency compared to Monte Carlo 
simulations which do not rely on this approximation. Figure 3(a) shows an example of exciton 
dynamics from the analytic model compared to Monte Carlo simulations with increasing numbers 
of injected excitons for a 5 /im long (9,8) SWCNT with diffusivity of 44 cm 2 /s and a linear relax- 
ation rate T = 3.3 ns -1 (i.e., %d = 5.7 ns), revealing good agreement. The quantum efficiency (Eq. 
(4)) for variable To can be calculated by Eq. (9) and agrees well with Monte Carlo simulations 
(Figure 3(b)). For a typical SWCNT, the effective quantum efficiency is significantly lower than 
the intrinsic quantum efficiency due to end quenching. The effective quantum efficiency is fur- 
ther reduced due to the highly efficient exciton-exciton annihilation, even for low exciton injection 
numbers (Figure 3(c)). 




100 200 300 400 500 1E-3 0.01 0.1 1 10 100 1 10 100 1000 

Time (ps) x D (ns) Injected Exciton Number (N ) 

Figure 3: (a) The exciton survival probability (Eq. (9)) for 1, 10, and 100 injected excitons on a 5 
/im SWCNT (diffusivity= 44 cm 2 /s, linear decay rate T = 3.3 ns -1 ) (b) The normalized quantum 
efficiency as a function of the average diffusional time with T = 3.3 ns -1 for ten injected excitons. 
(c) The normalized quantum efficiency as a function of injected exciton number for To = 5.7 ns 
and r = 3.3 ns -1 (corresponding to a 5 /im SWCNT with a diffusivity of 44 cm 2 /s). 

The analytic analysis (Eq. (9)) recreates the order of magnitude difference in saturation density 
observed in photoluminescence studies between suspended and encapsulated samples 2 -" 2 ^ ( Figure 
4(a) where PL(No, T£>,r) = NoTIqe(No : To,r)). The parameters used for encapsulated samples are; 
L = 610 nm, Dx = 4 cm 2 /s, T = 50 ns -1 and for air-suspended samples, L = 6 /im, Dx = 44 
cm 2 / s, r = 3.3 ns -1 . For comparison purposes, we define the saturation point (N sat ) with a simple 
model where PL(Nq) oc 1 — exp(—No/N sat ). From Figure 4(a) it can be seen that, as experimentally 
observed, the injected exciton density at which photoluminescence saturation occurs in "pristine" 
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SWCNT samples is about an order of magnitude than in a bulk encapsulated sample. Any re- 
maining discrepancies in the saturation fluence difference^ 2 -'- is likely due to random SWCNT 
orientation relative to incident light polarization, nonuniform excitation intensity, random length 
distribution, and different optical absorptions of bulk samples. 

As discussed earlier, even if the PL action cross-section is well quantified, it is challenging 
to decouple optical absorption from quantum efficiency. With fluences sufficiently high to cause 
exciton-exciton annihilation, the analytic model can be used to extract the optical absorption cross 
section (<J a b) f° r SWCNTs since it is directly linked to injected exciton number through the flu- 
ence, Nq = N c x o a t, x Fluence, where N c is the number of carbon atoms. For the most constrained 
results, independent measurements of the diffusivity (e.g., from PL line width), length (e.g., imag- 
ing), and decay lifetime (e.g., femtosecond excitation correlation spectroscopy ^J* 3 - 2 -), should be 
performed. As an example, Figure 4(b) shows a fit of the analytic model to experimental data 
for a single 4.6 /im (9,8) air-suspended SWCNT [from Xiao et al.— ]; the extracted absorption 
cross section is o a b = 5x 10~ 17 cm 2 /atom, which is about a factor of two higher than previously 
reported for SWCNTs on substrates. 33 Inputs to this model are the intrinsic linear lifetime of 218 
ps (extracted from fitting the model to FEC data, see below) and a diffusivity of 44 cm 2 /s.— 

The analytic model allows direct calculation of the underlying decay dynamics responsible for 
the sharp saturation in photoluminescence. For example, for each initial exciton number shown 
in Figure 4(c), the inset shows the corresponding normalized population decay dynamics com- 
pared to purely linear decay Pp. Increasing the injected exciton number significantly speeds up the 
population decay due to EEA, even at low injected exciton numbers. 

This analytic model also shows good agreement with photoluminescence lifetime experiments 
and can be used to extract the intrinsic decay rate T, without the need for additional fit parameters 
provided by a phenomenological bi-exponential decay model.— An example of this is shown in 
Figure 4(d) where normalized femtosecond excitation correlation data [see Xiao et al.— for ex- 
perimental details] is fit with the analytic model calculated for two injected pulses separated in 
time [see Supporting Information]. The complete model including end quenching and EEA shows 
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Figure 4: (a) Calculated photoluminescence for a typical encapsulated (red dashed curve) and 
air-suspended SWCNT (black solid curve) with the saturation point marked (as defined in text), 
(b) Analytic model (red curve), fit to experimental results (black circles) for a 4.6 /im (9,8) air- 
suspended SWCNT. (c) [INSET] Exciton population decay dynamics as a function of injected 
exciton number with purely linear decay Pr for comparison. Each exciton decay curve corresponds 
to a point on the photoluminescence vs injected exciton number plot, (d) FEC data (black circles) 21 
fit with model based on single exponential decay (blue dashed curve) and complete analytic model 
(red solid curve) and (bottom) corresponding residues. 
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better agreement to the data compared to simple linear decay. In both cases, only one fit parame- 
ter is used, namely T. Simple linear decay does not include end quenching, which is particularly 
relevant in the data for short time delays (<10 ps), causing the linear relaxation process to appear 
faster than it is (18.5 ns _1 for Pp compared to 4.6 ns _1 for Pjot)- We anticipate that the model 
proposed here will be particularly useful to the effort in understanding the length dependence of 
exciton lifetimes 14 since it avoids additional fit parameters. Comparison of the models shows that 
when end quenching is ignored the intrinsic exciton lifetime, and consequently the intrinsic dif- 
fusion length, is significantly underestimated. Knowledge of a sample's intrinsic properties will 
allow for engineering of optoelectronic properties such as length or diffusivity to maximize the 
quantum efficiency (for optical emission applications) or decay rate (for applications requiring fast 
relaxation) . 

Conclusion 

In this work an analytic model of exciton decay dynamics and quantum efficiency is developed 
including the three dominant exciton decay channels, linear decay, end/defect quenching, and 
exciton-exciton annihilation. Results from this analytic model are verified with Monte Carlo simu- 
lations which require considerable processing time for convergence. The spatial dependence of the 
exciton lifetime due to diffusion assisted quenching is explicitly included removing the need for 
approximations and phenomenological parameters. Consequently this model correctly accounts 
for the length dependence of the exciton lifetime and can therefore be used with SWCNT length 
studies to extract intrinsic diffusion lengths or with exciton lifetime studies to extract intrinsic 
exciton decay rates. 

In the multiple exciton regime, proper modeling of the exciton pair distribution function re- 
veals that for realistic exciton densities, the commonly used one dimensional ring approximation 
significantly underestimates EEA. Modeling including EEA and end effects reveals a method of 
extracting the optical absorption cross section for SWCNTs using photoluminescence saturation 
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experiments where the injected exciton number is directly linked to the optical pump fluence. From 
this it is seen that occurrence of photoluminescence saturation at an order of magnitude different 
excitation fluence, as observed in past experiments with samples with different diffusional times 
and linear lifetimes, is reasonable when both end effects and EEA are accounted for. 

Methods 

The Monte Carlo simulations are performed by considering a line segment of length L. Nq excitons 
are created at random positions along the interval [0,L] at time t = 0. For each time step St each 
exciton takes a step of length Sx = y/lDxSt in a random direction. Typically a value of St = 50 fs 
is used, however for small Xd a smaller St, down below 1 fs, is used to achieve convergence. The 
probability that an exciton decays via a linear pathway for each temporal step is TSt << 1. If an 
exciton leaves the segment [0,L] it is destroyed. If two excitons spatially overlap, one is destroyed. 
Generally the results were averaged over > 50,000 runs and each run continued until no excitons 
were left in the system. Numerical integration for the analytic model was performed using the 
Gauss-Kronrod quadrature method. 
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